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Abstract 

We report recent simulation results which might indicate the existence of a new low- 
temperature "phase" in an Ising lattice gas, driven into a non-equilibrium steady 
state by an external field. It appears that this "phase", characterized by multiple- 
strip configurations, is selected when square systems are used to approach the ther- 
modynamic limit. We propose a quantitative criterion for the existence of such a 
"phase" . If confirmed, its observation may resolve a long-standing controversy over 
the critical properties of the driven Ising lattice gas. 



1 Introduction 



Nearly two decades ago, Katz, Lebowitz and Spohn introduced a seemingly 
trivial generalization of the Ising model: a lattice gas driven far from equilib- 
rium by an external "electric" field, E [1]. The motivations are twofold: the 
statistical mechanics of non-equilibrium systems and the physics of fast ionic 
conductors [2] . Since then, this model has provided a rich variety of surprises, 
showing how our intuition built on the foundations of equilibrium statistical 
mechanics fails when applied to non-equilibrium stationary states. Though 
some of its unusual properties are now understood, it continues to offer new 
surprises [3]. In this paper, we report recent simulation results which may be 
signals of a novel "phase" . If this new state is established, it would lead to a 
resolution of a long standing controversy concerning the critical properties. In 
the next section, we briefly describe the model and the disagreement between 
Marro, et. al. [8] and the Leung- Wang team[10]. Since the methods used in 
these two approaches are quite distinct, we explore, in Section III, the pos- 
sibility that the discrepancies may lie in the existence of different types of 
ordering. In particular we introduce a new parameter which distinguishes the 
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two low temperature "phases" . The systems used in our simulations are sim- 
ple generalizations of the standard model - to include anisotropic interparticle 
interactions - so that the presence of the novel phase is more pronounced even 
for small lattices. The final section is devoted to some concluding remarks. 



2 The driven lattice gas and its critical properties 



The driven lattice gas is based on the Ising model [4,5] with attractive nearest- 
neighbor interactions. In the spin language, it is a ferromagnetic Ising model 
with biased spin-exchange [6] dynamics. Since both the Onsagcr solution [7] 
and most simulations concern a system in two dimensions {d = 2), we give a 
brief description of that model only. Each of the x Ly sites of a square lattice 
(on a torus, i.e., with periodic boundary conditions: PBC) may be occupied by 
a particle or left vacant. A configuration of our system is specified by the occu- 
pation numbers {rij}, where i is a site label and n is either 1 or 0. The interpar- 
ticle attraction is modeled by the usual Hamiltonian: Ti, = —4:JJ2<i,j>'ni'nj, 
where < i,j > are nearest-neighbor sites and J > 0. In thermal equilib- 
rium, a half filled system undergoes a second order phase transition, in the 
thermodynamic hmit, at the Onsager temperature Tq — (2.2692..) J//cb. To 
perform Monte Carlo simulations for this system, particles are allowed to hop 
to vacant nearest neighbor sites with probability min[l, ^-^'^/''bT-^^ where A7i 
is the change in Ti after the particle-hole exchange. The deceptively simple 
modification introduced in [1] is to bias the hops along, say, the y axis, so that 
the new rates are min[l, Q-i^'H+EAy)/kBT-^_ Locally, the effect of E is identical 
to that due to gravity. However, due to the PBC, this modification cannot 
be accommodated by a (single-valued) Hamiltonian. Instead, the system set- 
tles into a non- equilibrium steady state with a non-vanishing global particle 
current. For temperatures above some finite critical Tc, the particle density in 
this state is homogeneous, while below the system displays phase segrega- 
tion. Superficially, the driven system behaves like the equilibrium Ising model. 
With deeper probing, dramatic differences surface. For example, there is no 
question that its critical properties fall outside the Ising universality class. For 
recent reviews of this and other differences, we refer the interested reader to 
[3]. 

Focusing on critical properties here, let us review a long standing discrepancy 
between two sets of Monte Carlo results. While one set [10] shows data entirely 
consistent with field theoretic renormalization group predictions [11], another 
group [8] finds very different results, for which no viable theoretical analysis 
exists [12]. Technically, there are two serious differences between the simulation 
methods. In this section, we conjecture that different types of ordering may be 
responsible for the discrepancies. In view of the presence of "shape dependent" 
thermodynamics in the disordered phase of this system [13], we believe that 
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the scenario pictured below is quite plausible. Let us begin with a brief review 
of the controversy. 

Consider the (untruncated!) two point correlation function (in spin language): 

G{x,y) = {[2n{x,y) — 1] [2n(0, 0) — 1]), and its Fourier transform: S{k,p) = 
Y,x.yG{x,y)e^'^^^''^^^''^Py^'^y\ the structure factor. Due to the conserved dy- 
namics, S{0, 0) = 0, so that non-trivial large distance properties are found 
in, e.g., 5'(1,0) and S'(0, 1). In equilibrium systems with square geometries 
(Lj. — Ly — L), isotropy dictates that 5'(1,0) = 5'(0, 1), which approaches 
the static susceptibility in the thermodynamic limit for T > Tq- As T — Tq, 
this quantity diverges (as L'^^'' for finite systems at To). Below Tq, ordering 
is displayed as phase segregation, into a strip (of, say, a particle rich region) 
aligned with either the x- (or y-) axis, so that S'(0, 1) (or 5'(1,0)) L^. 
For the driven system, simulations [1] first revealed a discontinuity singular- 
ity (S'(1,0) > 5'(0, 1)) in the disordered phase and showed further that, at 
low T, phase segregation occurs only into "vertical" strips (S'(1,0) = O(L^), 
S'(0,l) = O(l)). Since these structure factors are related to (the inverse of) dif- 
fusion co-efficients in a continuum theory, these observations formed the basis 
for the study of a model where, first, the diffusion co-efficients are anisotropic, 
and second, more importantly, only one of these co-efficients vanishes as T ^ 
Tc [11]. This approach is also consistent with subsequent simulations [14] show- 
ing that, while S'(1,0) diverges as T ^ T,.. S{0, 1) remains 0(1) throughout 
the transition. A serious consequence of this "strongly anisotropic" behavior is 
that parallel and transverse lengths and momenta scale with different powers 
{Py ~ Px^^j A = 1 in mean field theory) in the critical region. Analogues of 
this unusual property exist in equilibrium systems [15]. When fluctuations are 
taken into account via renormalization group methods [11], the upper critical 
dimension (d^) is found to be 5, while anomalous scaling modifies A to 2 in 
d — 2. A more subtle result is that one set of exponents retains the classical 
values. 

Now, in order to /measure/ critical properties through simulations, with some 
accuracy, it is crucial to exploit finite size scaling methods [16]. Systems with 
various L's are used, and the quality of data collapse serves as an indicator 
of the presence of scaling. Given the highly anisotropic scaling of lengths, the 
analysis would be much simplified if samples obeying Ly ~ are used. This 
route is pioneered and followed in one of the studies [10]. More importantly, 
the order parameter used is closely related to 5'(1,0), so that it is sensitive 
to the presence of a single domain (i.e., one vertical strip, regardless of the 
system size). Also measured are cumulant ratios associated with S'(1,0). All 
data thus obtained display good collapse, using a set of exponents that is 
entirely consistent with those predicted by the renormalization group [11]. In 
another set of studies [8] , only square samples are used. Furthermore, the order 
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parameter used is 



(1) 



where M^^^^ 



Note that, if the ordered phase 



L ^ J2y{x) T^xiy) {"^IT-i " 1} 

consists of a single vertical strip, then m provides information about its average 
density. However, configurations with vertically aligned multiple strips also 
contribute to m. The conclusion is that, e.g., /3 ~ |, contradicting the field 
theory results [8]. 



3 A possible novel phase 



It should be quite alarming that two sets of simulations lead to drastically 
different conclusions, unless we allow the possibility that the "ordered" states 
in the two studies are not the same. Considering m is not sensitive to the 
presence of multiple strips, we conjecture the following. 

• In the thermodynamic limit, the low temperature state in systems with 
square geometry {L^ = Ly) is characterized by multiple strips, with a non- 
trivial distributions of widths. 

In other words, we believe that there is no long range order (which would lead 
to a single-strip state) while strip-strip correlations are controlled by a finite 
correlation length. For short, we will use a more picturesque term: "stringy 
phase." 




Fig. 1. A typical "stringy" configuration vs. a typical single strip one. (See text for 
system parameters.) 

Significantly, the existence and behavior of such configurations were first noted 
in [8]. Indeed, these authors wrote: "The system was never seen to escape 
from (these multiple-strip) states... and we had to manipulate some system 
configurations for L = 50 to create artificial one-strip states...". For L > 
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100 systems, they were "unable to observe the decay (of these multiple-strip 
states) toward one-strip states..." The conclusion was that multiple states are 
"metastable" with exceedingly long life times. 

Another source for our conjecture comes from simulation studies of a driven 
lattice gas with anisotropic interparticle interactions (J^ 7^ Jy) [17]. Though 
typical configurations are disordered/ordered (single strip) for high/low tem- 
peratures, there is a significant range of T where the system appears to be 
"stringy." In Fig.l, we show such a configuration at T = 0.80, found in a 
90 X 90 lattice with 3Jx = = J, driven with saturation E. (Here, as in 
the following, T is given in units of To.) In contrast, the single strip configu- 
ration was obtained with T = 0.70. To be more precise, the "stringy" states 
may be characterized by two properties: (a) the densities in each column (i.e., 
along the drive) are bi-modally distributed, yet (b) S{1, 0) remains small [18]. 
An example is provided in Fig. 2. While the column densities begin to develop 
a double peak around T ~ 0.90, the structure factor displays an inflection 
point at T ~ 0.70 (which is also where the fluctuations in S'(1,0) reach a 
maximum). Note that bi-modal column densities correspond to m > 0. Yet, 
in this range (0.75 ^ T < 0.95), the typical configurations are far from being 
single-stripped. To show that multiple-strip states are not simply metastable, 
we ran simulations starting from a completely ordered (single-strip) state. Un- 
like in the previous study [8], we do observe the split-up into multiple strip 
state^. In Fig.3, we show a time trace of 3500 measurements (700K MCS) 
of S{k,0), k = 1,2,3, which are sensitive to the presence of /c-strip domains. 
Note that S{1, 0) is often smaller than 5'(2, 0) and sometimes even smaller 
than S'(3,0), as shown in the inset. As a check, we made three separate runs 
of this length and observed the same type of behavior. Our conclusion is that, 
despite being started in a completely ordered state, the system evolves to- 
wards a "stringy" state. In other words, the probability to find the system in 
multiple strip configurations can be considerable. 

In order to make quantitative statements about the "stringy state" , we intro- 
duce the ratio: 



for a system in d dimensions. The numerator is a measure of "long range order" 
along the y-axis, while the denominator detects the presence of a single strip. 
In particular, we are interested in the behavior of 7?-, as L — 00, when the 
system is in various states. 

^ Figures in the Marro studies (e.g., Fig 2.7 in [9]) show typical configurations of 
T well below the transition value. It is possible that the failure to observe break- 
up/re- merging is due to T being too low (for the sizes studied). 
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(a) Histogram for column densities at T = 0.80[-], 0.90[«], 2.0[+], and (b) 
vs. T in a 30 X 30 lattice with 3Jx = Jy/'i = J-, driven with saturation E. 
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Fig. 3. Time trace of S{k^Q), k = 1 [solid line], 2 [•], 3, [dashed line], for a 60 x 60 
lattice at r = 0.85, driven with saturation E. Inset: magnified view of a portion. 

Far in the disordered phase, the two-point correlation decays as r~'^ [19]. Thus, 
G(0,L/2) 0{L-'^), while S{1,0) remains 0(1). On the other hand, deep in 
the ordered phase G 0(1) while S 0{L'^). As a result, we expect 



7^^0(l) 



(3) 



for both high and low temperatures. In between, a finite system makes a tran- 
sition, either through a "stringy phase" or directly into the ordinary ordered 
state (as in the equilibrium case). In the latter scenario, the results of the 
renormalization group analysis [11] should hold: G(0,r) ^~ld~2+A]/{i+A) ^^^^ 
S{k, 0) k~^. Thus, for simulations in = 2, we conclude that TZ decreases 
with L : 



7^ ^ 



(4) 
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In sharp contrast, for a stringy state, ordering has set in along y (so that 
G'(0, L/2) — 0(1)) but complete phase segregation in x is yet to take hold (so 
that 5(1, 0) — >• 0(1)). The consequence is an increasing TZ : 

n L'^. (5) 



which is a function of L. The drastic differences between these asymptotic 

properties (3, and 4 vs. 5) should be detectable through simulations. As a 
final remark, we note that, for the equilibrium Ising model, TZ decreases expo- 
nentially as L^e-^/S for T > T^, while for T <T^,n is 0(1). 

Turning to Monte Carlo data, we find that the evidence in favor of stringy 
states is quite strong, at least for the L's we can access. In particular, using 
models with J^, = J/3 and Jy = 3 J under saturation drive, we collected data 
from a number of systems with L e [10, 90]. The runs are up to 800K MCS 
long, with both random and ordered initial conditions. In Fig.4, we show a 
plot of TZ vs. L^, for various temperatures. For T = 0.85 and 0.90, this ratio 
appears to be proportional to L^. Meanwhile, its behavior is consistent with 
TZ —>■ 0(1) at both higher and lower T's. We believe that these results support 
the presence of a "stringy phase" . 
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Fig. 4. Correlation- Structure Factor ratio TZ for systems with a = 3 and saturation 
drive, plotted against L^. The legend refers to the temperatures of the various runs. 

Of course, this "phase" may be just a mirage in the landscape of a giant 
crossover. In the absence of analytic results, we can never be certain of its 
existence, since all simulations are based on finite L. Nevertheless, if future 
simulations continue to support 7Z oc 0{L'^), then it behooves us find a sound 
analytic basis for such a state. To conclude this section, we conjecture further. 
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As L increases further, the stringy phase will persist at lower and lower tem- 
peratures, so that it is the only low temperature state, for square systems, in 
the limit L — > oo. However, if the thermodynamic limit is approached through 
systems obeying Ly oc L^, the low temperature phase will consist of the or- 
dinary, single- domain state. If this conjecture turns out to be true, then an 
avenue is opened for the resolution of the controversy over critical properties. 



4 Concluding remarks 

In this short paper, we explored the possibility that the low temperature phase 
of the driven lattice gas depends on the geometry of the system. For square 
systems, we conjecture that a new, "stringy" phase is the only ordered state in 
the limit L — > oo. On the other hand, if an Ising-like, ordered state with a single 
domain is desired, then the "unusual" thermodynamic limit, with rectangular 
systems obeying Ly oc L\, must be used. For finite square geometries, we have 
some Monte Carlo evidence for the presence of a "stringy" state: multiple 
strips of high/low density regions. The system appears to be "ordered" in 
the direction parallel to the drive, but the strips are of varying widths and 
multiphcities. Specifically, as T is lowered, the disordered, homogeneous phase 
gives way to this multi-strip state for a sizable range of T, before settling into 
the fully phase segregated, single domain state at lower temperatures. 

Many open questions, besides those mentioned above, naturally arise. More 
quantitative analyses can be carried out, in order to better characterize this 
"stringy" state. For example, we could measure the distribution of strip- 
widths, or compile histograms for both the structure factors and the corre- 
lations. Once a clearer picture emerges, assuming such a state still exists, we 
should attempt to develop a reliable field theory. Certainly, the order param- 
eter of such a theory, m in Eqn. (1), is necessarily quite different from the one 
for an equilibrium Ising lattice gas. Hopefully, renormalization group methods 
can be applied, a nontrivial fixed point can be obtained, and its associated 
critical properties can be computed. Until then, to compare field-theoretic 
predictions with simulation data from the Marro studies makes as much sense 
as "comparing apples with oranges." 
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